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' in future, when more precise data on p meson form factors is available. 
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1 Introduction 



The QCD sum rules method [1] is one of the most powerful tools in studying low energy 
hadron physics. In this method, physically measurable quantities of hadrons are connected 
with QCD parameters, where hadrons are represented by their interpolating quark currents. 
The main idea of this method is to calculate the correlator functions of the interpolating 
quark currents in the deep Eucledian region with the help of operator product expansion 
(OPE) which allows one to take into account both perturbative and nonperturbative con- 
tributions. Relevant physical quantities are determined by matching the correlator to its 
phenomenological representation. 

In the current literature, a new, widely discussed alternative to traditional sum rules, 
namely, QCD light cone sum rules (QLCSR) is a convenient tool for the study of exclusive 
processes. This method is based on OPE on the light cone, which is an expansion over 
twists of the operators, rather than dimensions, as is the case in the traditional QCD sum 
rules. Moreover, in this method, all the nonperturbative dynamics encoded in the light cone 
distribution amplitudes, determine the matrix elements of the nonlocal operators between 
the vacuum and and the hadronic states (more about this method and its applications can 
be found in [2, 3]). 

In the present work, we study the p meson form factors in QLCSR. It should be men- 
tioned here that the p meson form factors are calculated at intermediate momentum transfer 
by using the three-point QCD sum rules method in [4]. Recently, p meson form factors, 
including next-to-leading-order perturbation theory, are analyzed within the same frame- 
work [5]. The Q 2 = point lies outside the applicability region of the three-point QCD 
sum rules. Consequently, extrapolating these form factors to Q 2 = 0, in principle, one can 
determine static characteristics of the p meson, such as magnetic and quadrapole moments. 
The p meson magnetic moment has already been investigated in the framework of the tra- 
ditional three-point and QLCSR in [6] and [7], respectively. It should be noted here that, 
the QLCSR is successfully applied to a wide range of problems of the hadron physics, for 
example, magnetic moments of the octet and decuplet baryons are calculated in [8] and [9], 
and magnetic moment of the nucleon in [10], respectively. 

The paper is organized as follows. In section 2, we derive the sum rules for the p meson 
electromagnetic form factors. In section 3, our numerical results and a comparison of them 
with the results of various approaches existing in the literature is presented. 

2 Theoretical framework 

In order to determine p meson electromagnetic form factors we will use QLCSR method. 
For this purpose, we consider the following correlator function 




(1) 



where 



J v (x) = u{x)~f v d(x) , 



q=u,d 
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are the p meson interpolating and electromagnetic currents, respectively. 

Firstly, we calculate the phenomenological part of the correlator. By inserting the 
complete set of states which have the quantum numbers of the p meson, between the 
currents in Eq. (1), we get 
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where p' = p + q, q is the momentum of the electromagnetic current, e is the polarization 
of the p meson and • • • describes contributions from higher states. The matrix element 
(p|j^|0) is determined as 



(p\ju\0) = f P m p el{p) 



(3) 



Assuming parity and time-reversal invariance, the electromagnetic vertex of the p meson 
can be written in terms of three Lorentz invariant form factors [11] 

(p(p,e)\J e J\p(P,e')) = e* a e'?{ - G^g^p + p% - G 2 (Q 2 )(q p9fia - q a g^) 

+ ^G 3 (Q 2 )q aqf3 (p + p')^ , (4) 

where Q 2 = —q 2 is the square of the momentum transfer. It should be mentioned here that, 
in practical computation, instead of calculating the Lorentz invariant form factors Gi(Q 2 ), 
the physical charge Gc, the magnetic Gm and quadrapole Gq form factors are often used. 

The Lorentz invariant form factors Gi(Q 2 ) are related to the charge, magnetic and 
quadrapole form factors through the following relations 
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Gn 
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G 2 , 
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where 77 = Q 2 /4m 2 . At zero momentum transfer, these form factors are proportional to the 
usual static quantities of charge e, magnetic moment p and quadrapole moment V 



eG c (0) = e , 
eG M (0) = 2m p p , 
eG Q (0) = m 2 V . 



(6) 



Substituting Eqs. (3) and (4) into Eq, (2), the phenomenological part of the correlator 
takes the following form 
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It is apparent from Eq. (7) that there are many different structures each of which can be 
used to extract the above-mentioned form factors. 

In order to find out the p meson electromagnetic form factors, we pick the following 
three structures: e u p p , £ nQ.v and (eq)p^q u . Hence we can write 



IV (p, q) = rii(p, q)e v p p + U 2 {p, q)e fJl q u + U 3 (p, q){sq)p^q u + ■■■ , (8) 



where 



Il 1 (p,q) = -G 1 (Q 2 )- 



2f p m p 



m 2 — (p + q) 2 
Ti 2 (p,q) = -G 2 (Q 2 )(l + 2r 1 ) 



m 2 - (p + q) 2 ' 

U 3 (p,q) = [2Gi(g 2 ) - G 2 (Q 2 ) + G 3 (Q 2 )(1 + 2 ^)] m 2 [m /!^ + qj] ' ^ 

For the invariant amplitudes Hi(p, q), one can write a general dispersion relation in (p + q) 2 
in the form 

TU(p,q)= J d8 a ^f q)2 +Ba bbr. , (10) 

where the spectral density corresponding to Eq. (7) is 

p l (s,q) = AF l (Q 2 )5(s-m 2 p ) , (11) 

where 

^2 



F 1 = -2G 1 (Q 2 )f p m p , 
F 2 = -G 2 (Q 2 )(l + 2r 1 )f p m p , 

F 3 = [2Gi(g 2 ) - G 2 (Q 2 ) + G 3 (Q 2 )(1 + 2 V )] 



fp m 



m 2 



According to QCD sum rules philosophy, in constructing sum rules for the form factors 
we need representation of the correlator function from QCD side. After contracting quark 
fields, the correlator function takes the form: 

IV = / rf 4 xe^(p(p)|e u M( a ;)7 M ^(x,0)7^(0) + e d M(0)7^ u (0,a;)7^(a:)|0) , (12) 

where S q (x) is the full quark operator of u (or d) quark. Imposing SU(2) symmetry and 
neglecting u and d quark masses, we have S u (x,0) = Sd(x,0) = S(x), and hence the full 
light quark propagator takes the form 



Six) = — — —r -(qq) — iq s / dv ^-t;G a g(vx) 7rr:VX a G a g'y3 
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Note that the second and third terms of S(x) do not give any contribution to the 
considered problem after Borel transformation is carried out (see below). Rewriting Eq. 
(13) in momentum representation and substituting it in Eq. (12), we get 

= j J rf 4 xrf 4 A;e^' ? -^|e tt (p(p)|M( a :)r i rf(0)|0)Tr7 M ^7,r i 
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where Tj = (J, 75, 7 Q , 7^75, a a p) is the full set of the Dirac matrices. It follows from this 
expression that, in order to calculate the correlator from QCD side one needs to know 
the matrix element of the nonlocal operators between vacuum and the p meson, i.e., p 
meson distribution amplitudes (or wave functions). We can see from Eq. (14) that main 
contribution to the correlator comes from only from the wave functions that contain odd- 
number of 7-matrices that are defined in the following way [14]: 
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where P^, e p and g pu are defined as 



where px = Px has been used. 

In Eqs. (15)— (18), (j)\\(v,/j,) is the leading twist-2 wave function, while g±, g\ and V are 
the twist-3 and all the remaining ones are the twist-4 wave functions. In Eqs. (17) and 
(18) the following relation is used 



where 



V(v,px) = J Vae ipx(ai+va:i) V(a u a 2 ,a 3 ) , 



T>a = daida2<ias5(l — ol\ — a 2 — a^) . 



In order to suppress the contributions of higher states and continuum and further elim- 
inate the subtraction terms in the dispersion relation, it is necessary to perform Borel 
transformation with respect to — (p + q) 2 . The contributions of higher states and contin- 
uum are subtracted using the quark hadron duality (for details, see [2, 3], [8]-[10], [12] and 
[13])- 

The sum rules for form factors can be obtained after applying Borel transformation 
to the two different representation of the invariant functions F t and then matching these 
results, and in doing so, we get 
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where M 2 is the Borel parameter and 



s(v) 



(Q 2 + m 2 v)v 



v — 1 — V , 
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2m 2 
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where Sq is the continuum threshold. 

In Eqs. (19)— (21), functions A\ and C u are defined as follows 



A i 



duA{u) , 

&( v ) = - r du$(u) , 

JO 
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duC(u) , 
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Note that the terms involving three-particle p meson wave functions are not presented 
in Eqs. (19)— (21) since their expressions are rather lengthy, however, their contributions 
are taken into account in the numerical analysis, which constitute about 5% of the total 
result. 



3 Numerical analysis 

In this section we present our numerical calculations on charge, magnetic and quadrapole 
form factors. The main input parameters of the QLCSR in regard to the above-mentioned 
form factors are the p meson wave functions, whose explicit expressions are given in [12] 
and [14], and we use them in our analysis. 

Apart from the wave functions, the sum rules for the form factors depend on the value of 
the continuum threshold Sq and Borel parameter M 2 . In the present work, we calculate the 
form factors at three different values of the threshold, i.e., s = 1.8 GeV 2 , s = 2.0 GeV 2 
and s = 2.2 GeV 2 (see also [4]). 

In Figs. (1), (2) and (3) we present the dependence of Gc, Gm and Gv on M 2 at 
six different values Q 2 = 0.5 GeV 2 , Q 2 = 1.0 GeV 2 , Q 2 = 2.0 GeV 2 , Q 2 = 3.0 GeV 2 , 
Q 2 = 4.0 GeV 2 and Q 2 = 5.0 GeV 2 of the momentum transfer, respectively. The Borel 
parameter M 2 in the sum rule is an auxiliary parameter and therefore physical quantities 
must be independent of it. For this reason we must determine the region of M 2 in which 
the form factors are independent of its value. In determining the working region of of the 
Borel parameter M 2 , the following two conditions must be satisfied: 

• M 2 should be large enough to suppress the contributions coming from higher twists, 
and, 

• M 2 should be small enough in order to suppress the continuum and higher state 
contributions. 

In the present analysis, both conditions are satisfied for the three form factors when M 2 
varies in the region 1.0 GeV 2 < M 2 < 2.5 GeV 2 . 
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In Figs. (4), (5) and (6) we present the dependence of the charge Gc, magnetic Gm and 
quadrapole Gq form factors on Q 2 , at three different values Sq = 1.8 GeV 2 , Sq = 2.0 GeV 2 
and s = 2.2 GeV 2 of the threshold, and at the fixed value M 2 = 1.0 GeV 2 of the Borel 
parameter. 

We observe from these figures that the dependence of the form factors on s is rather 
weak and when s vary from s = 1.8 GeV 2 to s = 2.2 GeV 2 . Unfortunately, the sum rules 
fail working at small Q 2 , and hence do not allow determination of the magnetic moment 
of the p meson with better accuracy compared to the results in the literature, such as, the 
results predicted by the sum rules [6, 7], Dyson-Schwinger based models [15, 16], Covariant 
light-front approach with constituent quark model [17], light-front formalism [18] and light- 
front quark model [19]. Our analysis predicts that, starting from Q 2 = 0.5 GeV 2 , the ratio 
Gm{Q 2 )/Gc{Q 2 ) is around 2.3 at all values of Q 2 . If we assume that this behavior holds at 
smaller values of Q 2 , we can conclude that p ~ 2.3, which is quite close to the prediction 
on the magnetic moment of the p meson [7]. It should be noted here that the values of 
the form factors at different values of Q 2 in this work are quite different compared to the 
predictions of the other approaches. For example, our results for Gc(Q 2 ) and Gm(Q 2 ) at 
all values of Q 2 are, approximately, two times smaller, while Gq(Q 2 ) is two times larger (in 
magnitude) compared to the prediction of [19] (we use the same parametrization for the 
form factors as in [19]). Therefore, more data on p meson is needed in order to choose the 
"right" model in calculating the form factors. 

In conclusion, we have presented the results for the p meson form factors in the frame 
work of QLCSR. We have obtained that, in the region of applicability of the method, at all 
values of Q 2 , the ratio Gm(Q 2 )/Gc(Q 2 ) is larger than 2, more precisely, this ratio varies 
around 2.3. This result can be checked when more precise data on p meson form factors is 
available. 
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Figure captions 

Fig. (1) The dependence of the charge form factor Gc on M 2 , at six fixed values of Q 2 , 
and at s = 2.2 GeV 2 . 

Fig. (2) The same as in Fig. (1), but for the magnetic form factor Gm- 

Fig. (3) The same as in Fig. (1), but for the quadrapole form factor Gq. 

Fig. (4) The dependence of the charge form factor Gc on Q 2 , at Three fixed values 
of s , and at M 2 = 1.0 GeV 2 . 

Fig. (5) The same as in Fig. (4), but for the magnetic form factor Gm- 
Fig. (6) The same as in Fig. (4), but for the quadrapole form factor Gq. 
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